Tensor force induced isospin-dependence of short-range nucleon-nucleon correlation 
and high-density behavior of nuclear symmetry energy 
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Quantitative information on the tensor force induced isospin-dependence of short-range nucleon- 
nucleon correlation (SRC) extracted from recent J-Lab experiments was used in constraining the 
high-momentum tails of single nucleon momentum distributions in both symmetric nuclear matter 
(SNM) and pure neutron matter (PNM). Its effects on the Equations of State of SNM and PNM 
as well as the nuclear symmetry energy are investigated. It is found that the tensor force induced 
isospin-dependence of SRC softens significantly the nuclear symmetry energy especially at supra- 
saturation densities. 
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I. INTRODUCTION 



Recently, much efforts have been devoted to study the 
isospin-dependence of short-range nucleon-nucleon cor- 
relations (SRC) in nuclei both experimentally 
theoretically [5l-tl0l|. For reviews, see, e.g., Refs. 
It is particularly exciting to note that experiments at 
the Jefferson Lab (J-Lab) have shown that about 20% of 
nucleons in 12 C are correlated. Most interestingly, the 
strength of np SRC is about 20 times that of pp (nn) 
SRC 3 . Theoretical studies have revealed that the dom- 
inance of np over pp (nn) SRC is a direct consequence 
of the tensor force acting in the np deuteron-like state 
@,0-0| • Moreover, it was found that [1, [l__] the isospin- 
dependence of SRC is robust and does not depend on 
the exact parameterization of the nucleon-nucleon force, 
the type of nucleus, or the exact ground-state wave func- 
tion used to describe the nucleons 0] . Because the local 
density of SRC pairs in nuclei is estimated to reach that 
expected in the core of neutron stars, it has been repeat- 
edly speculated that the observed isospin-dependence of 
the SRC may have significant effects on the Equation of 
State (EOS) of cold dense neutron-rich nucleonic matter 
and thus properties of neutron stars 0, 0. Neverthe- 
less, a quantitative evaluation of the effects has not been 
carried out yet. 

Nuclear symmetry energy E sym (p), which encodes the 
energy related to neutron-proton asymmetry in the nu- 
clear matter EOS, is a vital ingredient in the theoretical 
description of neutron stars and of the structure of neu- 
tron rich nuclei and reactions involving them. Since the 
density-dependence of E sym (p) is still the most uncer- 
tain part of the EOS of neutron-rich nucleonic matter 
especially at supra-saturation densities, to better deter- 
mine the E sym (p) has become a major goal of both nu- 

The E sym (p) can 



clear physics and astrophysics 
be approximated as the difference between the energy 
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per nucleon in pure neutron matter (PNM) and symmet- 
ric nuclear matter (SNM), i.e., E sym (p) = Epnm(p) — 
Esnm(p)- The isospin-dependence of SRC affects differ- 
ently the EOSs of PNM and SNM at high densities, it 
is thus expected to play an important role in determin- 
ing the high density behavior of E sym (p). In fact, it has 
long been known that the tensor force acting in isosinglct 
neutron-proton pairs influences the high density behav- 
ior of E sym (p) |27H3l| . However, the strength of the in- 
medium tensor force and its effects on the SRC were not 
so clear. The newly available and more quantitative in- 
formation on the tensor force induced isospin-dependence 
of SRC may allow us to better understand effects of the 
tensor force on the EOS and symmetry energy of dense 
neutron-rich nucleonic matter. In this work, within a 
phenomenological model for nuclear matter we show that 
the tensor force induced isospin-dependence of SRC stiff- 
ens significantly the EOS of SNM especially at super- 
saturation densities. However, it has almost no effect on 
the EOS of PNM. Consequently, the tensor force acting 
in neutron-proton isosinglet pairs significantly lowers the 
symmetry energy of dense neutron-rich nucleonic matter. 



II. NUCLEON MOMENTUM DISTRIBUTION 
INCLUDING EFFECTS OF SHORT-RANGE 
NUCLEON-NUCLEON CORRELATIONS 

To treat properly SRC effects has been a challenging 
task for nuclear many-body theories. For the present 
study, the most relevant question is how the SRC affects 
the one-body nucleon momentum distribution n(k). In 
fact, this question has been studied both experimentally 
and theoretically for a long time. Fortunately, many in- 
teresting results have already been well established and 
can be used reliably in our present work. For a com- 
prehensive review, we refer the reader to the book by 
Antonov et al. [32j . The n(k) contains information not 
only about the nuclear mean- field but also the SRC. For 
a Fermi gas of uncorrelated nucleons at zero temperature, 
the n(k) is simply a step function of one for k < kp, and 
zero for k > kp where kp is the nucleon Fermi momcn- 
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turn. The SRC, however, will deplete states below the 
Fermi surface and make the states above it partially oc- 
cupied, leading to a high-momentum tail in n(k). This 
feature has been well established in various microscopic 
and phenomenological models [32| . Thus, in a nucleus, 
the momentum of a low- momentum nucleon with k < kp 
is balanced by the rest of the nucleus; however, a high- 
momentum nucleon with k > kp is almost always bal- 
anced by only one other nucleon and the two nucleons 
form a correlated pair pj]. It is also well known that the 
SRC can be caused by either the repulsive core of the 
central force or the tensor part of the nucleon-nucleon 
interaction. The experimental finding that the np SRC 
dominates over the nn (pp) one indicates that the tensor 
force instead of the repulsive core is mainly responsible 
for the high-momentum tail of n(k). 
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FIG. 1: (Color online) Single nucleon momentum distribution 
in symmetric nuclear matter and pure neutron matter. 



To investigate how the strength of tensor force affects 
the high- momentum tail of n(k), some phenomenolo gica l 
methods have been shown to be particularly useful (32| . 
Instead of seeking for exact many-body wave functions 
containing the SRC as in microscopic many-body theo- 
ries, the phenomenological methods include the SRC in 
a physically very transparent way and allow one to eas- 
ily examine effects of the tensor force on the momentum 
distribution n(k). For instance, Dellagiacoma et. al have 
derived formulas for n(k) explicitly including the tensor 
force induced SRC for finite nuclei [H, H3 . It was shown 
that the probability of finding nucleons with high mo- 
mentum increases with the strength of tensor force (or 
cquivalently the percentage of D-wave mixture) [3^, H3] ■ 
With the same D-wave percentage, it was also found that 
the high-momentum tails of n(k ) for different nuclei are 
very close to each other [33|, Hj] . This finding is consis- 
tent with the conclusion of microscopic calculations that 
the behavior of n(k) at high momenta (k > 2 fm -1 ) is 
almost independent of the mass number [3^, [H| ■ More- 
over, the high momentum behavior of n(k) for symmetric 
nuclear matter was shown to be very similar to those of 
finite nuclei [H, Hf| HH • Results obtained using the y- 
scaling analysis of inclusive electron scattering data for 



2 H, 3 - 4 He, 12 C, 56 Fc also support this conclusion [13 • 

Here we parameterize the n(k) of nucleons in both 
SNM and PNM as 



n(k) = a (k < kp) 

= e bk {k>k F ), 
under the normalization condition (3o| that 
3 



F 



n(k)k 2 dk = 1. 
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(2) 
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The a and b are parameters determined using the exper- 
imental findings from the J-Lab experiments. For SNM, 
we take 80% (20%) nucleons as uncorrelated (correlated) 
as indicated by the J-Lab experiments and model calcu- 
lations as we discussed earlier. Thus, for SNM we have 
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n(k)k 2 dk = 0.8, 
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n(k)k 2 dk 



0.2. 
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The required parameter a = 0.8 and b is obtained nu- 



k 2 dk 



0.2. 



merically from the integration p- f£ 
The value of b depends on the Fermi momentum and 
thus the density. At nuclear matter saturation density, 
b = —1.899. Shown in the left panel of Fig[T]is the n{k) 
of symmetric nuclear matter at saturation density. The 
blue dashed line denotes the n(k) in an ideal Fermi gas. 
The red solid line is our parameterization. The black 
dotted line is the parameterization given by Ciofi degli 
Atti et.al. jsB], which is a fit to the result of variational 
many-body calculations [Hj]. It is clearly seen that our 
parameterization is very close to the microscopic single 
nucleon momentum distribution. For PNM, the SRC is 
induced only by the repulsive core. Utilizing the experi- 
mental finding that only 2% of nucleons in 12 C can form 
the nn- or pp- type SRC pairs, we have for PNM 
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—pr / n(k)k 2 dk = 0.98, 
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0.02. 
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The resulting parameter a is 0.98. Again, the parameter 
b is density-dependent and has to be found numerically. 
At nuclear matter saturation density, b = —2.435. The 
corresponding n(k) for PNM is shown in the right panel 
of Fig[T] As mentioned above, the big difference in n(k) 
for SNM and PNM is caused by the tensor-force induced 
SRC in SNM. 



III. EFFECTS OF SHORT-RANGE 
NUCLEON-NUCLEON CORRELATIONS ON 
THE EOS AND SYMMETRY ENERGY OF 
DENSE NEUTRON-RICH NUCLEONIC MATTER 

We now examine effects of the isospin-dependence of 
SRC on the EOS and symmetry energy of dense neutron- 
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rich nucleonic matter by using the single nucleon mo- 
mentum distributions in SNM and PNM constrained by 
the J-Lab experiments within a phenomenological ap- 
proach. With momentum-dependent nuclear interac- 
tions, not only the kinetic but also the potential en- 
ergy are affected by the high-momentum tail due to the 
SRC. Let us first examine separately effects of the SRC 
on the kinetic and potential energy per nucleon in both 
SNM and PNM. While the average kinetic energy per 
nucleon in a Fermi gas of independent nucleons is simply 

i~2m~' with correlated nucleons it is given by 
h 2 k 2 

— — n(k)k 2 dk, (5) 
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and for SNM and PNM, respectively. 



We compare in Fig[2] the average kinetic energies for 
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FIG. 2: (Color online) Average kinetic energy per nucleon 
Ekin for both symmetric nuclear matter and pure neutron 
matter. 

both SNM and PNM calculated with and without the 
SRC. As one expects, the SRC increases the Eu-n sig- 
nificantly for SNM but only little for PNM because of 
the strong isospin-dependence of the SRC. More quan- 
titatively, for SNM at saturation density, the E^m with 
the SRC (E kin {k° F ) ~ 40 MeV) is about twice of that 
(Ekin(kp) ~ 22 MeV) for the Fermi gas of uncorre- 
lated nucleons. However, for PNM, the increase in Ekin 
due to the SRC induced only by the repulsive core is 
very small. The different changes in Ekin for SNM and 
PNM due to the isospin-dependence of the SRC has a 
dramatic effect on the kinetic part of the nuclear sym- 
metry energy. In the Fermi gas without the SRC, the 
kinetic energy contribution to the symmetry energy is 

kin _ V Mn Ef%M = (2? - !)(§ ), which is B&- 



rpkin rpkin 

^sym — ^PNM 



ways positive with increasing density. However, because 
the SRC induced increase in Ekin is much larger for SNM 
than PNM, the kinetic contribution to the symmetry en- 
ergy now becomes negative. Thus, the kinetic part of 
E sym and its density-dependence are strongly affected 
by the tensor- force induced SRC. 
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FIG. 3: (Color online) Average potential energy per nucleon 
Epot for both symmetric nuclear matter and pure neutron 
matter. 



The average potential energy per nucleon E pot can 
be calculated using potential energy density function- 
als. For extracting information about the EOS of 
dense nuclear matter from heavy-ion reactions, several 
momentum-dependent single-particle potentials, such 
as the GBD (Gale-Bertsch-Das Gupta) HI, BGBD 
(Bombaci-Gale-Beretsch-Das Gupta) |39j and MDI 
(Momentum-Dependent Interaction) [4(J phenomenolog- 
ical potentials, have been widely used. In calculating the 
corresponding EOSs with these potentials only the step 
function for the momentum distribution of a Fermi gas 
was used. Thus, effects of the SRC were not considered. 
Here we use the MDI interaction and the nucleon momen- 
tum distributions including the SRC effects. The MDI 
potential energy per nucleon can be written as (40| 



Epot — — 



2 Po 
C 
PPo 



D 



P 



lPo 



(0) 



n(ki)n(k 2 ) 



1 + (fci - k 2 ) 2 /A 2 



- dk\ dki , 



where a = 4/3 and the parameter A = l.Ofc^ [40j. The 
n(ki) and n{k,2) are the one-body momentum distribu- 
tion of nucleon-1 and nucleon-2. For SNM, as usual, the 
parameters A = -114.079 MeV, B = 107.154 MeV and 
C = —0.0274 MeV are obtained from fitting its three em- 
pirical saturation properties, namely the vanishing pres- 
sure, the average binding energy of E — —16 MeV, and 
the incomprcssibility of Kq = 220 MeV at the saturation 
density po = 0.16/m~ 3 . For PNM, we obtain the three 
parameters A = -120.570 MeV, B = 102.080 MeV and 
C = —7.61 x 10~ 4 MeV by satisfying the well established 
theoretical constrains on the low density PNM EOS given 
in Refs.[4l|,|42[, namely, E = 1.7 MeV at k F = 0.4 fin" 1 
and E = 4.2 MeV at k F = 0.8 fm , as well as the sym- 
metry energy E sym = 31 MeV at p . Shown in Fig[3] 
is a comparison of E pot with and without the SRC for 
both SNM and PNM. Similar to the case of E kin , the 
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SRC increases the magnitude of E pot significantly for 
SNM but has negligible effect for PNM. Consequently, 
the tensor force induced SRC also softens the potential 
part of the symmetry energy = E% M - Ef^ M . 

Shown in Figg] are the EOSs of SNM and PNM with 
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the isosinglet nucleon-nucleon interaction dominates the 
symmetry energy [43|, [44| . It was also known that the 
tensor force acting in the isosinglet channel affects sig- 
nificantly the high-density behavior of nuclear symmetry 
energy [27l - [3]| . The finding here that the SRC softens 
the symmetry energy is thus not surprising qualitatively. 
In fact, the short-range repulsive tensor force may even 
lead to negative symmetry energies at supra-saturation 
densities as first predicted by Pandharipande and Garde 
in 1972 within a variational many-body approach (27| . 
However, we never had before the kind of quantitative 
information on the isospin-dependence of the SRC due 
to the tensor force and the repulsive core from experi- 
ments. The results presented here thus represent a sig- 
nificant progress in our understanding about the tensor 
force effects on the EOS and symmetry energy of dense 
neutron-rich nucleonic matter. 



IV. SUMMARY 



FIG. 4: (Color online) Total binding energy per nucleon for 
both symmetric nuclear matter and pure neutron matter. 
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FIG. 5: (Color online) Symmetry energy with (correlated) 
and without (uncorrelated) the shore-range nucleon-nucleon 
correlations. 

and without the SRC. The corresponding symmetry en- 
ergies are presented in FigjSJ It is interesting to see that 
the tensor force induced SRC has a significant impact 
on the symmetry energy especially at supra-saturation 
densities. The symmetry energy obtained in calculations 
with the SRC is much softer especially for kp > 2.0 fm -1 . 
Microscopic many-body theories have shown clearly that 



In summary, the quantitative information on the ten- 
sor force induced isospin-dependence of SRC from the 
recent J-Lab experiments was used in constraining the 
high-momentum tails of single nucleon momentum dis- 
tributions in SNM and PNM. The latter are then used 
in evaluating the EOSs of SNM and PNM as well as the 
symmetry energy of neutron-rich nucleonic matter. It is 
found that the isospin-dependence of SRC stiffens signif- 
icantly the EOS of SNM but has little effect on that of 
PNM. Consequently, the nuclear symmetry energy espe- 
cially at supra-saturation densities is significantly soft- 
ened. This finding confirms the critical role of the tensor 
force acting in neutron-proton isosinglet pairs in deter- 
mining the high-density behavior of nuclear symmetry 
energy predicted by various microscopic and phenomeno- 
logical nuclear many-body theories. 
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